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Abstract 

Multi-physics  simulations  based  on  multi-component  multi-solver  modeling  approach  were  performed  for  high-temperature  fuel  cells.  The 
developed  approach  was  primarily  aimed  at  the  design  of  complex  multi-component  engineering  systems.  It  extends  the  libraries  of  ear¬ 
lier  designed  multi-physics  system  with  component  classes,  which  makes  it  particularly  suitable  for  modeling  of  fuel  cell  systems.  The  C++ 
based  class  hierarchy  enables  simple  implementations  of  different  physical  models  based  on  general  3D  PDE  (partial  differential  equations) 
solvers,  or  simplified  engineering  ID  or  2D  models.  Simulations  of  solid-oxide  fuel  cells  were  performed  using  a  combined  transport  solver 
in  multi-species  environment.  The  components  included  the  PEN  complex  (anode,  cathode,  electrolyte),  air/fuel  channels,  interconnects,  seals 
and  ambient  environments.  Species  concentrations,  mass,  momentum,  energy  fluxes,  and  electric  potentials  were  solved  for  different  compo¬ 
nents.  Models  for  unsteady  fluid  dynamics  of  species,  heat  transport,  electrochemistry  and  electric  currents  were  combined  within  different 
components  and  interfaced  for  common  variables  at  the  inter-component  boundaries.  The  results  include  steady  and  unsteady  distributions 
of  temperature,  species  concentrations,  mass  fluxes  and  electric  potential  inside  co-flow  and  cross-flow  fuel  cells  with  different  number  of 
channels. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Background 

Computer  simulations  of  fuel  cells  involve  complex  multi¬ 
physics  modeling,  and  earlier  work  on  modeling  these  systems 
was  faced  with  the  challenge  of  unifying  different  solution  pro¬ 
cedures  [1-3].  The  complexity  of  modeling  is  due  to  the  fact 
that  fuel  cells  are  both  multi-physics  and  multi-component  sys¬ 
tems.  The  components  of  a  typical  fuel  cell  are:  anode,  cathode, 
electrolyte  (PEN  complex)  as  well  as  interconnect,  seal,  and 
air/fuel  channels.  Different  physical  processes  occur  inside  each 
component:  electrochemistry  (PEN),  mass  transport  (channels), 
electrical  current  (PEN,  interconnects),  heat  transfer  and  associ¬ 
ated  thermal  expansion,  etc.  In  a  continuous  effort  of  design 
optimization  these  components  are  arranged  in  a  variety  of 
geometrical  designs,  such  as  tubular,  planar,  hexagonal,  and  in 
various  air/fuel  flow  configurations:  co-flow,  counter  flow,  cross¬ 
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flow,  as  well  as  mixed  configurations.  In  addition  to  this,  fuel 
cells  arrangements  in  stacks  also  present  a  complex  optimization 
problem. 

In  the  previous  work  the  authors  developed  a  multi-physics 
approach  for  modeling  a  broad  class  of  problems  of  discrete  and 
continuum  dynamics  [4],  which  was  implemented  in  a  multi- 
physics  simulation  system MulPhys  (http://www.mulphys.com). 
This  approach  is  especially  suitable  for  handling  segregated  or 
multi- component  systems,  where  different  physical  phenomena 
occur  in  different  regions  of  space.  A  variety  of  naturally  occur¬ 
ring  and  human-made  structures  fall  into  this  category,  such  as 
biological  and  geological  structures,  as  well  as  many  engineering 
systems.  In  this  study  the  approach  was  applied  to  the  modeling 
of  fuel  cells,  which  represent  typical  multi-component  multi- 
physics  systems.  In  modeling  these  systems  it  is  important  to 
provide  a  simulation  environment  capable  of  incorporating  dif¬ 
ferent  solution  methods  and  numerical  schemes,  while  enabling 
their  inter-dependence.  A  distinct  feature  of  a  fuel  cell  system 
is  that  some  physical  processes,  such  as  gas  flow  and  electrical 
fields,  are  restricted  to  certain  regions  of  space,  while  other,  such 
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as  heat  transfer  span  multiple  regions.  Our  challenge  in  design¬ 
ing  the  simulation  system  was  to  incorporate  this  feature  into 
the  modeling  paradigm. 

2.  Method 

Earlier  simulations  of  fuel  cells  performed  by  the  authors 
[5-7]  prompted  the  introduction  of  a  physical  component  as 
a  new  entity  of  the  multi-physics  modeling  framework  [4].  A 
component  is  this  case  represents  a  well  defined  region  of  space 
governed  by  one  or  more  physical  laws.  Fig.  1  shows  the  UML 
(Unified  Modeling  Language)  diagram  [8]  (http://www.uml.org) 
of  the  class  hierarchy  designed  for  the  adequate  representation  of 
multi-component  systems.  In  this  description  the  system  consists 
of  domains ,  where  each  domain  is  a  collection  of  components , 
with  each  component  represented  by  a  set  of  models.  Each  model 
is  defined  by  a  collection  of  variables  and  methods  to  manipulate 
them. 

In  a  single  processor  simulation  one  domain  is  usually  used  to 
represent  one  fuel  cell,  while  in  a  multi-processor  case  multiple 
domains  are  used,  with  each  domain  assigned  to  one  proces¬ 


sor.  In  this  case  one  domain  can  still  represent  a  single  fuel 
cell  or  a  part  of  it.  Multi-processor  simulations  can  be  used 
to  model  the  behavior  of  stack  of  fuel  cells,  or  to  do  a  statis¬ 
tical  analysis  of  their  performance,  such  as  risk  analysis.  In 
the  first  case  a  fast  inter-processor  communication  is  essen¬ 
tial  for  the  efficiency  of  simulation,  and  a  computer  cluster 
is  usually  an  appropriate  implementation  medium.  In  the  sec¬ 
ond  case  a  statistically  independent  sample  is  needed,  and  no 
inter-processor  communication  is  necessary.  A  large  number  of 
distributed  loosely  coupled  processors,  such  as  provided  by  an 
emerging  grid-computing  environment  will  be  suitable  for  this 
purpose. 

Following  this  approach,  the  setup  of  a  system  consists  of 
two  parts:  geometric  design  and  physical  modeling.  During 
geometric  design  the  shape  of  the  domain  is  created,  includ¬ 
ing  all  its  components.  Then  each  component  is  discretized 
into  control- volumes  or  elements.  Each  element  will  be  auto¬ 
matically  assigned  a  number  corresponding  to  the  type  of  the 
component  it  belongs  to. 

In  the  stage  of  physical  modeling  each  model  is  defined  by 
specifying  dependent  variables  and  methods  to  update  them 


Fig.  1.  UML  diagram  of  class  hierarchy. 
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Fig.  2.  Design  schematics  for  a  multi-physics  multi-component  system. 


(iterators).  Models  represent  physical  processes  and  are  not 
related  to  any  particular  geometry  or  part  of  a  physical  space.  In 
the  final  stage  the  components  are  specified,  with  each  compo¬ 
nent  being  a  simple  list  of  models  which  belong  to  it.  At  this  point 
the  geometry  is  linked  to  physics  by  assigning  a  set  of  models  for 
each  component  (Fig.  2).  This  scheme  enables  an  easy  incorpo¬ 
ration  of  complex  modeling  of  multi-physics  multi-component 
systems. 

Each  model  is  a  member  of  a  model  class  and  is  character¬ 
ized  by  a  set  of  variables,  and  two  functions:  the  constructor,  or 
initializer  and  the  iterator .  Within  each  model  there  is  an  access 
list  to  the  elements  which  belong  to  the  model  can  be  seen  in  Fig. 
2,  this  linking  elements  to  models  is  done  indirectly  through  the 
components.  At  the  same  time  there  is  a  direct  link  from  each 
model  to  domain  elements.  This  scheme  allows  one  to  assign 
multiple  models  to  each  element,  thus  enabling  multiple  physics 
in  the  same  region  of  space.  At  the  same  time,  component-based 
modeling  also  enables  one  to  restrict  some  physical  processes 
to  specific  regions  of  space.  For  example,  if  the  heat  transfer 
is  occurring  in  multiple  components  of  the  domain,  the  heat 
transfer  model  will  have  a  list  of  pointers  to  all  the  elements 
of  the  domain  involved  in  heat  transfer.  On  the  other  hand,  the 
gas  transport  model  will  be  present  only  in  components  associ¬ 
ated  with  air/fuel  channels.  The  element- to-component  links  are 
generated  during  the  geometric  design  of  the  system,  while  the 
model-to-element  lists  are  automatically  created  during  the  call 
to  model  initializer  function. 

Both  initializer  and  iterator  can  access  various  geometrical 
information  for  each  element.  This  information  is  either  directly 
stored  for  each  element,  or  is  computed  on-the-fly.  For  exam¬ 
ple,  the  dimensionality  of  the  element,  it’s  order,  number  of 
faces,  and  the  access  list  to  the  vertexes  of  the  element  is  directly 
stored  in  memory  for  each  element.  On  the  other  hand,  the  infor¬ 
mation  on  the  face-center  vectors,  face  areas,  and  face  normal 
vectors,  etc.  can  be  computed  from  the  stored  information.  In 
some  applications  it  may  be  favorable  to  change  the  ratio  of 
stored  to  computed  data  for  greater  efficiency.  For  example, 
the  face  areas  and  normal  vectors  can  be  computed  during  the 
initialization  stage  and  stored  in  memory  for  the  duration  of 
the  simulation.  At  the  same  time,  in  the  case  when  the  geom¬ 


etry  is  changing  with  time  these  values  should  be  computed 
on-the-fly. 

It  should  be  noted  that  a  general  control- volume  discretiza¬ 
tion  scheme  used  in  this  modeling  framework  can  be  used  for 
both  detailed  differential  as  well  as  coarse-grained  integral  anal¬ 
ysis.  A  differential  analysis  can  involve  finite-difference,  finite 
volume  or  finite  element  schemes  for  solving  PDEs,  while  an 
integral  analysis  can  be  based  on  ID  or  2D  engineering  approx¬ 
imations  and  simplified  control- volume  balance  schemes.  The 
latter  has  an  advantage  of  simplicity  and  efficiency.  For  example, 
the  fluid  flow  in  micro-channels,  such  as  those  used  in  fuel  cells, 
can  be  approximated  with  a  Poiseuille  flow  solution,  instead  of 
resorting  to  a  computationally  expensive  finite  volume  or  finite 
element  approximations.  Similarly,  a  2D  surface  electrochem¬ 
istry  model  can  be  used  instead  of  a  more  accurate,  but  more 
expensive  differential  solution  inside  the  PEN  complex.  Such 
an  approach  opens  a  possibility  for  multi-scale  multi-physics 
modeling  in  a  simplified  control-volume  framework. 

This  modeling  framework  was  implemented  in  C++  language 
on  a  Beowulf  cluster.  The  setup  and  control  of  the  simulations  on 
the  cluster  was  done  remotely  from  a  workstation  using  a  Java- 
based  GUI,  which  provides  the  features  of  remote  login  and 
monitoring  of  distributed  simulations,  and  capable  of  complex 
geometric  design  [9-11]. 

2.7.  Mathematical  formulation 

In  the  fuel  cell  modeling  scheme  four  basic  sub-models  are 
used:  mass  transport,  species  transport,  heat  transport,  and  elec¬ 
tric  current.  The  mass,  momentum,  energy  and  species  conser¬ 
vation  equations  are: 


7  [  pdV  +  [  p(V  ■  n)dA  =  -mwaii, 

ot  Jv  J A 

VpdV  +  J  Vp(V  -n)dA  =  XA 

—  J  ^totPdV  +  J  (e  +  —  +  — )  PV  w)dA  =  Qm, 

7T"  [  pYk  dV  +  [  pYk(V  -n)dA  =  cbkAwa h 
ot  J\/  J  A 


where  etot  is  the  total  stored  energy  per  unit  mass,  e  is  the  internal 
energy  per  unit  mass,  and  (9 in  is  the  net  heat  flux  into  the  volume. 
Wall  mass  and  charge  fluxes  are: 


fi^wall  —  FwaiiAwan/9,  CQfc  —  (2) 

nkF 

where  n  is  the  number  of  electrons  per  model  of  reactant  k.  The 
PEN  and  separator  plates  are  considered  to  be  made  of  solid 
material,  therefore  only  the  energy  equation  was  solved  in  these 
regions,  which  was  simplified  to: 


d 

dt 


^totpdV  —  Gin  ~b  Ggen 


where  the  radiative  and  convective  heat  fluxes  through  the  sur¬ 
face  of  the  control  volume,  and  the  thermal  energy  transported 
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by  mass-flux,  are  all  included  in  Gin?  and  the  heat  source,  Ggen, 
is  obtained  from  ohmic  heating.  The  heat  of  entropy  generation 
results  in  the  following  expression: 

•  V  ^  •  •  •  • 

Gin  =  /  v  Gcond  T  G  conv  +  Grad  +  G  mass? 

A 

Ggen  —  Oden)-^  —  TAsC0\{2 

The  total  entropy  change  per  mode,  As  is  obtained  from 

rp 

As  =  As0  +  Ru  In  — 

rP 

where  As0  is  the  molar  change  in  entropy  at  standard  conditions 
and  tr  and  rp  are  the  reactant  and  product  activities  respectively. 


2.2.  Numerical  discretization 

The  spatial  discretization  method  used  in  the  solver  is  based 
on  control- volume  techniques,  and  is  applicable  to  complex 
geometries.  The  basic  spatial  discretization  elements  are  of 
an  arbitrary  polyhedral  type,  which  makes  them  easily  adapt¬ 
able  to  arbitrary  shaped  domains,  and  generally  produces  better 
quality  meshes  than  those  based  on  simple  tetrahedral  or  hexa- 
hedral  types.  The  geometry  input  for  the  solver  was  generated 
using  a  voxel-based  graphical  user  interface  developed  by  the 
authors  [9-11],  and  specifically  designed  for  describing  multi- 
component  structures. 

The  discretized  equation  for  the  total  internal  energy  can  be 
obtained  from  the  energy  equation  (1): 


Ainode  (  O2  ^cathode  \  ^ 

rR  -  —  (  ^  J 


rP 


^H20  ^anode 


Pressure,  P  is  calculated  from  the  ideal  gas  law: 

P  =  pRuT 

The  electrochemistry  model  is  based  on  the  assumption  that 
the  overall  chemical  reaction  occurring  in  the  fuel  cell  is: 

H2(g)  +  502(g)  =>■  H20(g) 


Calculation  of  the  cell  potential  starts  with  the  Nernst  equa¬ 
tion: 


E  =  E°  +  —  (  In 
IF 


[Xh2][Xq2]1/2 

[Xlbo] 


The  pressure  is  assumed  to  be  the  same  for  both  the  anode  and 
cathode  gas  channels.  The  reversible  potential  at  standard  state 
conditions  is  obtained  from  the  change  in  the  standard  Gibbs 
free  energy. 


AG° 

nF 


The  corrected  cell  potential,  Ecov,  is  obtained  by  subtract¬ 
ing  the  ohmic  ( Oohm ),  concentration  (1 7con),  and  activation  (r]act) 
losses  (i.e.  overpotentials)  from  the  ideal  Nernst  potential,  E : 


Vohm  7?  con  0  act 


All  overpotentials  are  related  to  the  current  density,  and  the 
activation  over-potential  is  defined  by  an  empirical  relation  rep¬ 
resented  by  a  limiting  form  of  the  Butler- Volmer  equation: 


7? ohm  —  Gen^net?  Ocon  — 


RUT  /  _  'den 


nE 


il 


^u  T  ,  (  Gen 

^7  act  =  - In 


anF 


G 


The  convective  heat  transfer  coefficient,  hc ,  is  calculated  from 
the  Nusselt  number,  which  is  a  function  of  the  Reynolds  and 
Prandtl  numbers:  hc  =  -. 


en  ■  on 

n+ 1  _  tot  ,in 

Got,  i 


Nf 


Pi 


n 


n.  V 

^  1  k=  1 


P/V/ 


(4) 


where,  the  summation  is  done  over  all  the  faces,  Nf,  of  a  poly¬ 
hedral  cell.  After  neglecting  potential  energy,  the  total  internal 
energy  can  be  found  from  the  absolute  enthalpy: 

p  V2  V2 

Got  j  =  G  4 - b  ~zr  ~  hi  +  ~zr 

Pi  2  2 

The  absolute  enthalpy  can  be  found  from 


h(T)  =  h°\Ttef  +  j  CpdT, 

Tcf 


\n 


Q'm,i  ~  2cond,(  +  Gconv.j 
Nf 


_l_  nn  1  nn  1  nn 

i  '  >Zwall ,i  '  >Zmass ,/  '  *Zs,i  ’ 


2 kki  Aj(T  —  Ti) 


k  +  ki  Ax, 


Gcond,  i  =  E 
i= 0 

Nf 

Gconv.i  =  E^^tot.il m  <  0  +  mkem,k\mi  >  0), 
(=0 

Nf 

Gwall.i  =  E  hMTk  -  T),  Qn 


_  •  n  n  /r\ 

mass,/  m  wall,  /Got,/ 


/=0 


Mass  fraction  for  species  s  =  1 , . . . ,  Ns  —  1  is  found  from 
species  transport  equation: 


yW+ 1 
1  s,i 


KiPNi  a? 


X 


p’Ni  Pi 

-  Nf 

•  nk)  +  qnDs,ipAk 

lk=  1 


d>s  ( Ms  A Wall,/ 


_ 

^ s,i 


—  l 


n 


den,/ 


zF 


(6) 


tfbs,i,k  ~  PDs 


yn  _ yn 

1  s,k  1  s,i 

Ax 


(7) 
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the  mass  fraction  for  species  Ns  can  be  found  from  the  definition 
of  mass  fraction: 


N*-l 


Yn+j  =  \~T  Y". 

ivs,t  /  j  s,i 


n+ 1 


s=  1 

The  molecular  weight  for  the  mixture  is  then  obtained  from: 

1 


M"+1  = 


Zsh(Y:f)/(Ms) 

The  molar  concentration  of  species,  X,  can  be  found  from: 


vU  1  _  y m + 1  i _ 

sj  ~  sj  Ms 

The  partial  pressures  for  each  species  is  then  obtained  from 
the  mole  fraction  and  the  total  pressure: 

v«+ 1 
pn-\~  1  _  _ 

$•>  i  pn  -f- 1 

ri 

The  mixture  density  is  obtained  from  the  ideal  gas  law: 


p',+1  = 


p 


n- |- 1 


i 


pt i+l  y  w+l 


where  the  gas  constant  for  the  mixture  can  be  calculated  as: 


R"+l  = 


Ru 


M 


n+ 1 


The  exit  velocity  can  be  found  from  the  mass  conservation 
equation: 

1 


vnJrl .  = 
exit,; 


At 

exit  ,iPi 


n+ 1 


X 


n+lw.  _  n/iw. 

Vi  vz  «  /f/W  -  \  4  .  « 

^  Pk  \yk  ’  nk)  Ak  ^wall,; 


yvf 


At 


&=i 


7Vf 


mwall 


,z  —  ^s,  i  ^ wall,  i 


(8) 


s=\ 


The  pressure  drop  is  obtained  from  the  momentum  equation: 


A  P 
Ax 


Nf 


V i~P> ?Vi  V-  n  n  -n  - 
'  '  -^kPm  ■  nk)Ak 


At 


k=i 


+fp«)2 


(9) 


where  the  friction  coefficient,  Cf  is  the  function  of  the  Reynolds 
number:  Cf  =  f(Re),  Re  =  pvd/ti. 


3.  Results  and  discussion 


The  most  important  processes  to  be  included  in  fuel  cell  mod¬ 
eling  are  heat  transfer  and  species  transport.  These  models  were 
implemented  on  the  basis  of  the  discretized  energy  transport 
equation  (4)  and  species  transport  equation  (6).  The  heat  trans¬ 
fer  scheme  was  validated  on  the  analytical  solutions  for  constant 
gradient  heat  transfer,  and  the  conservation  of  mass  and  momen¬ 
tum  has  been  guaranteed  through  (8)  and  (9). 


Fig.  3.  Geometry  of  a  cross-flow  fuel  cell. 


A  planar  fuel  cell  was  investigated  for  two  typical  flow  con¬ 
figurations:  co-flow  and  cross-flow.  The  geometry  of  a  planar 
fuel  cell  consists  of  a  set  of  straight  channels  separated  by 
a  PEN  complex  (anode-electrolyte-cathode)  and  surrounded 
by  seals  and  interconnects.  The  5  air  x  5  fuel  channel  con¬ 
figuration  was  investigated  for  both  co-  and  cross-flow  cases. 
Fig.  3  shows  a  typical  cross-flow  geometry  with  5x5  chan¬ 
nels,  which  was  produced  using  a  specialized  voxel-based  GUI 
[9-11]. 

A  general  control-volume  discretization  employed  in  this 
modeling  framework  was  used  for  a  coarse-grained  integral 
analysis  of  flow,  heat,  and  species  transport  in  various  fuel-cell 
geometries.  Fig.  4  shows  two  cross-sections  of  the  computa¬ 
tional  grid  of  13  x  13  x7  control  volumes  for  a  5  x  5  co-flow 
fuel  cell.  Different  components  are  shown  by  different  colors, 
which  include  air/fuel  channels,  PEN,  interconnects,  and  chan¬ 
nel  inlet/outlet  regions. 

Figs.  5  and  6  show  the  temperature  distribution  in  air  and  fuel 
channels  of  a  co-flow  fuel  cell.  The  temperature  distribution  in 
the  PEN  of  a  co-flow  fuel  cell  is  shown  in  Fig.  7.  The  temperature 
increase  in  the  direction  of  the  flow  is  caused  by  accumulated 
heat  release  in  the  PEN  complex  and  is  typical  for  a  co-flow  con¬ 
figuration.  These  figures  should  be  compared  to  the  respective 
temperature  distributions  for  the  PEN  (Fig.  8)  and  air/fuel  chan¬ 
nels  obtained  for  the  cross-flow  fuel  cell  (Figs.  9  and  10).  The 
uneven  distribution  is  due  to  the  localized  heat  release,  which 
occurs  in  the  areas  where  both  air  and  fuel  channels  have  a  con¬ 
tact  through  a  PEN  component.  The  temperature  distribution  is 
also  asymmetric  with  respect  to  either  air  or  fuel  channel  axes, 
which  is  typical  for  the  cross-flow  geometry. 

Fig.  11  shows  species  distribution  along  the  air/fuel  chan¬ 
nels.  As  can  be  seen,  even  though  the  fuel  is  consumed,  oxygen 
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Fig.  4.  Cross-sections  through  a  5  x  5  cross-flow  fuel  cell. 


concentration  experiences  no  appreciable  change  because  of  a 
relatively  high  mass  flow  rate  in  the  air  channel. 

Fig.  12  shows  the  convergence  of  the  temperature  to  steady 
state  by  considering  the  variation  of  three  points  located  in  the 
center  of  the  fuel  cell  with  time.  The  three  points  are  located 
in  the  fuel  channel,  PEN,  and  air  channel  at  the  center  of  the 
five  cell  cross-flow  cell.  The  thermal  response  time  of  the  cell 
is  much  longer  then  the  response  of  other  phenomena  occurring 
in  the  cell  thus  it  is  used  to  identify  whether  the  cell  has  reached 
steady  state. 

Temperature  and  species  distributions  for  different  fuel  cell 
configurations  are  in  a  qualitative  agreement  with  the  behavior 
of  typical  fuel  cells.  In  particular,  in  a  previous  study,  Zitney  et  al. 
[12]  coupled  a  FLUENT  CFD  model  with  Aspen  Plus  (a  process 
simulator)  to  simulate  a  single  channel  planar  SOFC  fuel  cell 
auxiliary  power  unit.  Stevenson  et  al.  [13]  and  Cheng  et  al.  [14] 
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Fig.  5.  Temperature  (in  K)  in  air  channels  of  a  co-flow  fuel  cell. 
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Fig.  6.  Temperature  in  fuel  channels  of  a  co-flow  fuel  cell. 
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Fig.  7.  PEN  temperature  distribution  in  a  co-flow  cell. 
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Fig.  8.  PEN  temperature  in  a  cross-flow  cell. 


Fig.  11.  Species  distributions  in  air/fuel  channels. 
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Fig.  9.  Temperature  distribution  in  air  channels  of  a  cross-flow  cell. 


conducted  similar  studies  and  presented  temperature  distribu¬ 
tion  for  co-flow,  counter-flow,  and  cross-flow  planar  SOFC.  The 
results  of  the  last  two  studies  show  that  the  co-flow  predicted  the 
most  uniform  temperature  distribution  while  the  cross-flow  was 
the  most  non-uniform.  The  cross-flow  temperature  distributions 
in  all  studies  had  a  similar  high  temperature  region  located  in  the 
corner  opposite  the  air  and  fuel  inlets.  Exactly  the  same  trends 


Fig.  10.  Temperature  distribution  in  fuel  channels  of  a  cross-flow  cell. 


Time  [s] 

Fig.  12.  Transient  temperature  at  three  control  points. 


and  the  range  of  temperatures  were  observed  in  this  study.  The 
temperature  distribution  in  PEN  (Fig.  8)  also  agrees  well  with 
the  one  obtained  by  Achenbach  for  a  cross-flow  fuel  cell  [1]. 

The  simulations  were  performed  on  a  Beowulf  cluster  of 
1  GHz,  1  GB  computing  nodes,  with  each  node  of  the  cluster 
running  a  simulation  of  a  single  cell.  The  speed  of  execution 
was  limited  by  the  process  with  the  minimum  time  scale  (electro¬ 
chemistry).  For  the  5x5  channel  cross-flow  case  a  convergence 
to  the  steady  state  took  about  2000  s  of  physical  time,  with  a  typ¬ 
ical  time  step  of  0.01  s,  which  corresponds  to  about  30  min  of 
CPU  time. 

4.  Conclusions 

A  multi-component  approach  to  the  modeling  of  multi¬ 
physics  systems,  such  as  fuel  cells  has  been  proposed  and  tested 
on  several  cases  of  planar  high- temperature  fuel  cells.  Cou¬ 
pled  with  the  appropriate  geometrical  design  tool,  the  approach 
enables  easy  incorporation  of  various  physical  models  for  dif¬ 
ferent  regions  of  space,  as  well  as  their  co-existence  in  the 
same  region.  In  particular,  the  approach  enables  coupled  solu¬ 
tion  based  on  numerical  schemes  for  3D  PDEs,  as  well  as  on  ID 
and  2D  engineering  approximations.  Temperature  and  species 
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distributions  for  different  fuel  cell  configurations  are  in  a  good 
agreement  with  the  earlier  studies  and  is  consistent  with  the 
behavior  of  a  typical  fuel  cell. 

The  approach  is  general  and  can  be  applied  to  the  modeling 
of  other  multi-component  systems  common  in  science  and  engi¬ 
neering.  For  example,  biological  tissues  and  geological  layers 
are  typical  segregated  media  where  the  current  approach  can  be 
effectively  used  for  complex  multi-physics  modeling. 

Parallel  implementation  of  the  approach  for  modeling  of  fuel 
cell  stacks  is  straightforward  and  is  the  subject  of  the  future 
work. 
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